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ABSTRACT 

Extensive calculations of the atomic data required for the spectral modelling 
of the K-shell photoabsorption of oxygen ions have been carried out in a multi- 
code approach. The present level energies and wavelengths for the highly ionized 
species (electron occupancies 2 < < 4) are accurate to within 0.5 eV and 0.02 
A, respectively. For A^ > 4, lack of measurements, wide experimental scatter, 
and discrepancies among theoretical values are handicaps in reliable accuracy 
assessments. The radiative and Auger rates are expected to be accurate to 10% 
and 20%, respectively, except for transitions involving strongly mixed levels. Ra- 
diative and Auger dampings have been taken into account in the calculation of 
photoabsorption cross sections in the K-threshold region, leading to overlapping 
lorentzian shaped resonances of constant widths that cause edge smearing. The 
behavior of the improved opacities in this region has been studied with the xs- 
TAR modelling code using simple constant density slab models, and is displayed 
for a range of ionization parameters. 

Subject headings: atomic processes — atomic data — line formation — X-rays: 
general 



- 2 - 



1. Introduction 

The unprecedented spectral resolution of the grating spectrographs on the Chandra 
and XMM-Newton X-ray observatories have unveiled the diagnostic potential of oxygen K 
absorption. O VII and O VIII edges are observed in the spectra of many Seyfert 1 galaxies 
which have been used by Lee et al. (2001) to reveal the existence of a dusty warm absorber 
in the galaxy MCG -6-30-15. Steenbrugge et al. (2003) have detected inner-shell transitions 
of O lll-O VI in the spectrum of NGC 5548 that point to a warm absorber that spans three 
orders of magnitude in ionization parameter. Moreover, Behar et al. (2003) have stressed 
that, in the case of both Seyfert 1 and Seyfert 2 galaxies, a broad range of oxygen charge 
states are usually observed along the line of sight that must be fitted simultaneously, and 
may imply strong density gradients of 2-4 orders of magnitude over short distances. 

Schulz et al. (2002) have reported that the spectrum of the archetypical black hole 
candidate Cyg X-1 shows a prominent oxygen K edge that when modelled suggests a deficient 
O abundance approximately consistent with an abundance distribution in the interstellar 
medium (ISM). Observations of the O I edge and the equivalent widths of O Il-O VIII in the 
spectrum of the black-hole binary LMC X-3 by Page et al. (2003) lead to upper limits of 
the neutral and ionized column densities that rule out a line-driven stellar wind, and thus 
suggest a black body with a multi-temperature disk. A gravitationally redshifted O Vlll Lya 
observed in absorption in the spectrum of the bursting neutron star EXO 0748-676 shows, 
unlike Fe ions, multiple components consistent with a Zeeman splitting in a magnetic field 
of around 1-2x10^ G (Loeb 2003). Sako (2003) has shown that in optically thick conditions 
an O VIII-N VII Bowen fiuorescence mechanism, caused by the wavelength coincidence of the 
respective Lya (2p-ls) and Ly(^ (7p-ls) doublets, can suppress O emission while enhancing 
C and N thus leading to abundance misinterpretations. 

The ls-2p resonant absorption lines from O ions with electron occupancies iV < 3 are 
being used with comparable efficiency in the detection of the warm-hot intergalactic medium 
which, according to simulations of structure formation in the Universe, contains most of the 
baryons of the present epoch (Cagnoni et al. 2004). A spectral comparison of low and high 
extinction sources can lead to the separation of the ISM and instrumental components of 
the K edge (de Vries et al. 2003), the former closely resembling the edge structure of neutral 
oxygen computed by McLaughhn & Kirby (1998) with the i?-matrix method (Berrington et 
al. 1987; Seaton 1987). High-resolution spectroscopy of the interstellar O K edge in X-ray 
binaries has been carried out by Juett et al. (2004), discounting oxygen features from dust 
and molecular components and providing the first estimates of the O ionization fractions. 

Accurate laboratory wavelength measurements for K-shell resonance lines in O Vll have 
been reported by Engstrom & Litzen (1995) using laser pulses focused on compressed powder 
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targets and by Beiersdorfer et al. (2003) with an electron beam ion trap; the latter method 
has been recently extended to O vi and O V (Schmidt et al. 2004). K-shell photoabsorp- 
tion of O I has been measured by Menzel et al. (1996) and Stolte et al. (1997), and of 
O II by Kawatsura et al. (2002) where the role of theoretical data — namely level energies, 
wavelengths, 5f/-values, and natural widths — in experimental interpretation has been em- 
phasized. Laboratory energies for K-vacancy states in O I and O ll have been obtained by 
Auger electron spectrometry (Caldwell & Krause 1993; Krause 1994; Caldwell et al. 1994), 
but worrisome discrepancies with the photoabsorption measurements have not as yet been 
settled. 

Such experimental activity has certainly encouraged theoretical verifications which have 
been carried out with a variety of atomic structure and close-coupling methods. Radiative 
and Auger data related to the sateUite lines of H-, He-, and Li-like oxygen have been com- 
puted in an approach based on the XjZ hydrogenic expansion (Vainshtein & Safronova 1971, 
1978) and with the multiconfiguration Dirac-Fock method (Chen 1985, 1986). A theoretical 
study of the K-shell Auger spectrum of O I has been carried out with the multiconfigura- 
tion Hartree-Fock method (Saha 1994) finding good agreement with the measurements by 
Caldwell & Krause (1993). Energy positions and total and partial Auger decay rates for 
the [ls]2p^ "^P and terms of O ll have been determined by Petrini & Araujo (1994) with 
the structure code SUPERSTRUCTURE (Eissner et al. 1974) and photuc (Saraph 1987), the 
latter being based on an implementation of the close-coupling approximation of scattering 
theory. The K-shell photoionization cross section of the ground state of neutral oxygen 
has been computed with the i?- matrix method by McLaughlin & Kirby (1998), giving a 
detailed comparison with the experimental energy positions for the [ls]2p^(^P)np and 
[ls]2p"^(^P)np "^P" Rydberg series. More recently, the i?-matrix method has also been used 
to compute the high-energy photoionization cross section of: the ground state of O VI in an 
11-statc approximation (Charro et al. 2000); O I taking into account core relaxation effects 
and the Auger smearing of the K edge (Gorczyca & McLaughlin 2000); and in its relativistic 
Breit-Pauh mode, O l-O VI where wavelengths and /-values are listed for the n — 1 Ka 
resonance complexes with / > 0.1 (Pradhan et al. 2003). 

We are interested in a realistic modelling of the oxygen K edge and in formally estab- 
lishing its diagnostic potential. For this purpose we have computed as accurately as possible 
the required atomic data, some of which have not been previously available for the whole 
oxygen isonuclear sequence, incorporating them into the XSTAR photoionization code (Kall- 
man & Bautista 2001). Previous models have been fitted with two absorption edges and five 
Gaussian absorption lines (Juett et al. 2004). The present data sets have been computed in 
a systematic approach previously used for the iron isonuclear sequence (Bautista et al. 2003; 
Palmeri et al. 2003a,b; Bautista et al. 2004; Mendoza et al. 2004) which has allowed the 
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determination of the efficiency of Fe K line emission and absorption in pfiotoionized gases 
(Kallman et al. 2004). Extensive comparisons are performed with several approximations 
and previous data sets in order to size up contributing effects and to estimate accuracy 
ratings. With the improved opacities, the edge morphology dependency on the ionization 
parameter is studied for the first time. 



2. Numerical methods 

The numerical approach has been fully described in Bautista ct al. (2003). The atomic 
data, namely level energies, wavelengths, gf -values, radiative widths, total and partial Auger 
widths, and total and partial photoionization cross sections are computed with a portfolio of 
pubhcly available atomic physics codes: AUTOSTRUTURE (Badnell 1986, 1997), HFR (Cowan 
1981), and bprm (Berrington et al. 1987; Seaton 1987). Although relativistic corrections are 
expected to be negligible for the lowly ionized species of the O sequence, they can be more 
conspicuous for members with greater effective charge z — Z — N + 1, Z being the atomic 
number and N the electron occupancy. For consistency, wavefunctions are calculated unless 
otherwise stated with the relativistic Breit-Pauli Hamiltonian 

-f^bp = -f^nr + -f^lb + H2h (1) 

where iifnr is the usual non-relativistic Hamiltonian. The one-body relativistic operators 

N 

Hib = /n(mass) + /„(d) + /„(so) (2) 

n=l 

represent the spin-orbit interaction, /„(so), the non-fine-structure mass variation, /„(mass), 
and the one-body Darwin correction, /n(d). The two-body Breit operators are given by 

H2h = 5Z 5'nm (so) + g„m{ss) + S'nm (cSs) + g„m{d) + g„m{oo) (3) 

where the fine-structure terms arc g,i,n{so) (spin-othcr-orbit and mutual spin-orbit), gnm{ss) 
(spin-spin), and the non-fine-structurc counterparts gnm{css) (spin-spin contact), gnm{d) 
(two-body Darwin), and gnmioo) (orbit-orbit). 

BPRM is based on the close-coupling approximation whereby the wavefunctions for states 
of an N-electron target and a colliding electron with total angular momentum and parity 
Jtt are expanded in terms of the target eigenfunctions 

^''^-'^J^x.^ + E^A-- (4) 

i 3 
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The functions Xi vector coupled products of the target eigenfunctions and the angu- 
lar components of the incident-electron functions; Fi(r) are the radial part of the continuum 
wavefunctions that describe the motion of the scattered electron, and A is an antisymmetriza- 
tion operator. The functions are bound-type functions of the total system constructed 
with target orbitals. Breit-Pauli relativistic corrections have been implemented in BPRM by 
Scott & Burke (1980) and Scott & Taylor (1982), but the inclusion of the two-body terms 
(see Eq. 3) is at present in progress. Auger and radiative dampings are taken into account 
by means of an optical potential (Robicheaux et al. 1995; Gorczyca & Badnell 1996, 2000) 
where the resonance energy with respect to the threshold acquires an imaginary component. 
In the present work, the A'"-electron targets are represented with all the fine structure levels 
within the n = 2 complex, and level energies are set to those obtained in approximation 
AS2. Positions for states of the {N + l)-clectron system are obtained from the peaks in the 
photoionization cross sections. This data set is referred to as RMl. 

Core relaxation effects (CRE) are investigated with AUTOSTRUTURE by comparing an 
ion representation where all the electron configurations have a common basis of orthogonal 
orbitals, to be referred to hereafter as approximation ASl, with one where each configuration 
has its own set, approximation AS2. Data sets computed with hfr are labelled HFl. 

3. Energy levels and wavelengths 

We have looked into several interactions that can influence the accuracy of the atomic 
data, namely configuration interaction (CI), relativistic corrections, and CRE. It has been 
found by calculation that CI with the n = 3 complex gives rise to small contributions that 
can be generally neglected. Computations are carried out in intermediate coupling that take 
into account relativistic corrections even though the latter are small in particular for low 
ionization stages. Some effort is then focused in characterizing CRE. 

Energies calculated with approximations ASl, AS2, HFl, and RMl for levels within the 
n = 2 complex of O ions are listed in Table 1. Since the accuracy of the transition data, e.g. 
wavelengths and gf -values, depends on the representations of both the lower (valence) and 
upper (K-vacancy) levels, detailed comparisons are carried out for each type. For the valence 
levels (see Table 2), a comparison of ASl and AS2 reveals the importance of CRE which 
lower energies by as much as 20%. Since RMl also neglects CRE, the best accord is with 
ASl: within 15% although the RMl energies are on average lower by 4%. It may be noted 
that the RMl calculation for O I has been performed in LS coupling due to the large size 
of the relativistic option. AS2 and HFl agree to 5% and within 10% with the spectroscopic 
values (Moore 1998) which we find satisfactory. 
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Regarding K- vacancy levels, the energy differences of the approximations considered 
with respect to AS2 are plotted in Fig 1. The comparison with ASl shows that, for species 
with a half-filled L shell or with greater electron occupancies (6 < iV < 8 and E < 544 eV), 
CRE again lower energies but by amounts that increase with A^: AE = E{AS1) —E{AS2) < 
0.7 eV for = 6 to AE ^ 2.5 eV for A^ = 8. On the other hand, for ions with a half-empty 
L shell or with lower occupancies (2 < A^ < 5 and E > 544 eV), level energies are raised by 
CRE particularly for A" = 3 where AE > —1.35 eV. RMl behaves, as expected, in a similar 
fashion to ASl, but some inconsistencies (e.g. lowered rather than raised levels) occur for 
N < 6 perhaps due to resolution problems in determining the position of narrow resonances 
in the photoionization curves. The agreement between the AS2 and HFl energies is a sound 
±0.6 eV which sets a limit to the accuracy level that we can attain, and encourages us to 
adopt AS2 and HFl as our working approximations. 

AS2 and HFl K-vacancy level energies are compared with experiment and other theo- 
retical data in Table 3. It may be noted that there are no reported measurements for species 
with 4 < A^ < 6. HFl level energies are within 0.5 eV of the spectroscopic values for A^ < 3 
and the CI values by Charro et al. (2000) for N = 2 while AS2 is undesirably higher (0.9 
eV) for the [ls]2p '^P° level in O VII. A convincing comparison with experiment for A^ > 7 is 
impeded by the scatter displayed in measurements which is well outside the reported error 
bars, and the photoionization measurements by Stolte et al. (1997) are 1 eV lower than 
HFl. The energy positions computed by Pradhan et al. (2003) with bprm are noticeably 
higher than present results for A^ > 4, in particular for the B-hke system (A" = 5) where 
discrepancies as large as 5 eV are encountered. The values listed by McLaughhn & Kir by 
(1998) for O II (A^ = 7), which are also distinctively high (~3 eV), correspond to the CI 
target they prepared for the non-relativistic i?-matrix calculation of the photoionization of 
O I. They contrast with those quoted for an equivalent target by Gorczyca & McLaughlin 
(2000) which are in good agreement with Stoltc ct al. (1997). However, the value obtained 
by McLaughlin & Kirby (1998) for the [ls]2p^ hole state in O I is in good agreement 
with AS2 and HFl. 

CRE also have an impact on the transition wavelengths where, as depicted in Fig. 2, the 
comparison of ASl and RMl with AS2 indicates that in general CRE lead to increasingly 
longer wavelength for A^ > 6 (A > 23 A) and shorter for higher charge states. It may also be 
seen in Fig. 2 that wavelengths obtained with AS2 and HFl agree to ±0.02 A setting again 
an accuracy bound. (AS2 and HFl wavelengths are also listed in Table 4.) The present 
wavelengths are compared with experiment and other theoretical predictions in Table 5. For 
A^ < 4 the present data, in particular HFl, agree with measurements and Vainshtcin & 
Safronova (1978) to ±0.02 A; in this respect, for 2 < A^ < 4 the values by Chen (1985, 1986) 
are somewhat long while those by Behar & Kahn (2002) are short, and for A" = 3 the values 
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by Pradhan et al. (2003) are also long. For 5 < N < 6 the outcome is less certain due to 
absence of laboratory measurements and inconsistencies in the computed values, e.g. the 
order of the upper levels in iV = 5. For > 6 comparisons with measurements seem to 
indicate that all theoretical results are slightly short although the agreement of HFl is still 
within ~0.02 A except for the transition from the ground level of O II and the measurement 
by Stolte et al. (1997) in O I which are just out (less than O.OSA). It must be pointed out 
again the undesirable scatter of the experimental data for O I. 



4. gf -values 

In the following discussion of (^/-values, the transitions involving the strongly mixed 
[ls]2p"^ ^0,3/2 and ^P3/2 levels in O II are excluded due to extensive cancellation that result 
in completely unreliable data; also the Is^ ^Sq [ls]2p '^P'j' intercombination transition in 
O VII with a very small gf-vahie {gf < 10~^) that would require a much greater effort to 
guarantee respectful accuracy. 

A comparison of 5f/-values obtained with the ASl and AS2 approximations shows that 
CRE are relatively unimportant, leading to differences much less than 10% except for tran- 
sitions involving the four [ls]2s2p ^P° relatively mixed levels of O VI where they jump up to 
~15%. In Fig. 3 we compare the AS2 and HFl gf-yalues (also hsted in Table 4) showing 
an agreement within 10% except for O I where the latter are consistently 25% higher. In a 
comparison of AS2 and HFl /-values with other theoretical data sets in Table 6, it is found 
that those by Pradhan et al. (2003) are generally 20% smaller; however, an interesting point 
arises regarding the discrepancies between AS2 and HFl where Pradhan et al. (2003) favors 
AS2 for O I and HFl for the small /-value in O VI (see Table 6). 



5. Radiative and Auger widths 

The radiative width of the jih level is defined as 

Aj — ^2 ^ji (5) 

i 

where Ajj is the rate (partial width) for a downward (j > i) transition. It is found by 
calculation that CRE disturb the radiative widths of K-vacancy levels in the O isonuclear 
sequence by 10% or less (5% for N > 5) except when strong admixture occurs, namely for 
[ls]2p^ ^D3/2 and ^P3/2 in O ll and the four [ls]2s2p ^P° levels in O VI where differences are 
around 20%. 
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By contrast, as shown in Fig. 4, CRE reduce Auger rates by amounts that grow with N: 
from 10% for N = 4 where log^a < 14 up to 35% for = 8; however, this trend is broken 
for N = 3 due to strong level mixing where those for the [ls]2s2p('^P°) ^Pj are increased by a 
factor of 2 (not shown in Fig. 4) even though the [ls]2s2p(^P°) ^Pj levels are hardly modified 
(less than 5%). The Auger widths obtained with HFl for the former levels are ~20% higher 
that AS2 while the remaining, as shown in Fig. 4, are on average lower by 10%. With this 
outcome we estimate an accuracy rating for the Auger rates at just better than 20%. As 
further support for this assertion, in Table 7 we compare our branching ratios for KLL Auger 
transitions with other theoretical and experimental values. It may be seen that the ASl, 
AS2, and HFl data are stable to better than 10%, and the accord with other data sets is for 
most cases better than 20%. Radiative and Auger widths computed in approximations AS2 
and HFl are listed in Table 8. 



6. Photoionization cross sections 

High-energy photoionization cross sections of O parent ions with electron occupancies 
3 < A" < 8 have been computed in approximation RMl. Intermediate coupling has been 
used throughout except for O I because of computational size and negligible relativistic 
corrections. Following Gorczyca & Robicheaux (1999) and Gorczyca & McLaughlin (2000), 
damping is taken into account in detail in order to bring forth the correct K threshold 
behavior. In Fig. 5 our cross sections are compared with those computed by Pradhan et al. 
(2003) with the bprm package in the region near the n = 2 resonances and those by Reilman 
& Manson (1979) in a central field potential where it may be seen that, even though the 
background cross sections are in satisfactory accord, the positions and sharpness of the 
edges are discrepant. In the present approach, due to the astrophysical interest of both the 
n = 2 resonances and the edge region, the resonance structure is treated in a unified manner 
showing that a sharp K edge occurs only in O VI due to the absence of spectator Auger 
(KLL) channels of the type 

[ls]2p''np ^ 2p'*-2np + e- (6) 
^ [2s]2p^-^np + e~ (7) 
[2s]22p^np + e- . (8) 

For the remaining ions, these KLL processes dominate over participator Auger (KLn) decay 



[ls]2p^np ^ 2p^-^+e- 
^ [2s]2p'^ + e- 



(9) 
(10) 
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causing edge smearing by resonances with symmetric profiles of nearly constant width. Par- 
ticipator Auger decay, on the other hand, gives rise to the familiar Feschbach resonances — 
displayed by O VI in the region below the K threshold (see Fig. 5) — whose widths become 
narrower with n thus maintaining sharpness over the spectral head. 

In order to elucidate further the interesting properties of K resonances, we analyze the 
experimental and theoretical predictions reported for the two [ls]2p'^(^P)np and [ls]2p^(^P)np 
Rydberg series in O I that exhibit significant discrepancies. Regarding resonance positions, 
a stringent test is to compare quantum defects, /j,, defined as 

E = E,,^ - -^^^ (11) 

where E^^ is the series-hmit energy in Rydbergs; z = Z — N + 1 is the effective charge 
{z = 1 for O l); and n is the principal quantum number. In Fig. 6, the RMl quantum 
defects for these two series are plotted relative to their respective series limits and compared 
with those obtained in photoionization measurements (Stolte ct al. 1997) and the i?-matrix 
calculation of McLaughlin & Kirby (1998). Firstly, it may be seen that the experimental 
error bars grow sharply with n; hence greater accuracy for high n is obtained by quantum- 
defect extrapolation than actual measurement. Secondly, the agreement of experiment with 
RMl is very good while that with McLaughlin & Kirby (1998), who used practically the 
same numerical method (i?- matrix), is only a poor 40%. This outcome may be due to their 
neglect of damping effects or their choice of a Hartree-Fock basis for the Is, 2s, and 2p 
orbitals which results in a deficient representation for the K-vacancy states. Since the R- 
matrix method employs a common set of target orbitals for both valence and K-vacancy 
states, we have found the most adequate orbital basis to be that optimized variationally 
with an energy sum comprising all the states in the target model. The importance of K- 
vacancy state representations has also been stressed by Gorczyca & McLaughlin (2000). 
Thirdly, it is found that the RMl quantum defects (n > 3) are hardly modified by small 
target-energy shifts prior to Hamiltonian diagonalization, a standard i?-matrix procedure to 
enhance accuracy by using experimental thresholds; therefore, the absolute energy positions 
of these resonances mostly depend on the actual position of their series limits. Moreover, 
as shown in Fig. 7, the huge and incorrect quantum-defect decline obtained for high n in 
McLaughlin & Kirby (1998) when using two different sets of experimental thresholds, namely 
those of Krausc (1994) and Stolte ct al. (1997), arc not due in our opinion to the threshold- 
energy shifts but to resonance misidentifications where /x(n -|- 1) — fi{n) ^ —1 instead of 
approximately zero (see Fig. 7). The quantum-defect changes for n < 6 in McLaughlin 
& Kirby (1998) as thresholds are shifted are very small in accord with our findings. By 
contrast, experimental quantum defects are markedly sensitive to the series-limit positions; 
for instance, if the quantum defects in the photoionization experiment of Menzel et al. (1996) 
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are computed with the series hmits by Krause (1994), they display noticeable discrepancies 
(see Fig. 7) while if they are computed with those by Stolte et al. (1997) (not shown), 
consensus is reached. 

It is inferred from the above comparison that the energy positions obtained by Auger 
electron spectrometry (Krause 1994; Caldwell et al. 1994) for the [ls]2p^ and states 
of O II and listed in Table 3 should perhaps be scaled down by 0.4 eV. If this proposition is 
extended to the [ls]2p^ '^P° state of O I, then the measurements by Stoltc ct al. (1997) and 
Krause (1994) arc brought to close agreement; the only remaining experimental discrepancy 
for the latter state would then be the measurement by Menzel ct al. (1996) still high by 1 eV. 
Furthermore, contrary to the [ls]2p'^('^P)np '^P° resonances of O I in our RMl calculation, 
the absolute energy position of [ls]2p^ ^P° is insensitive to threshold-energy shifts. This is 
a consequence of its wavefunction being mostly represented by the second term of the close- 
coupling expansion (4). Therefore, a useful parameter in comparisons with experiment — in as 
much as it gives an indication of both a balanced close-coupling expansion and experimental 
consistency — is the energy interval between the lowest n = 2 and n — 3 resonances, AE{2, 3) 
say. In Table 9, the RMl AE{2, 3) for the whole O isonuclear sequence are tabulated and 
compared with experimental and theoretical estimates. The 0.07 Ryd (1 eV) discrepancy 
in the A£'(2, 3) obtained from the measurements by Stolte et al. (1997) and Menzel et al. 
(1996) in O I (A^ = 8) denotes irregular data: since the positions for the (^P)3p state differ 
by only 0.07 eV the issue is then the [ls]2p^ ^P°. The agreement of RMl with Stolte et 
al. (1997) is very good; however, if the RMl (^P) threshold is improved by replacement 
with the experimental value at Halmiltonian diagonalization, then AE{2, 3) is reduced to 
1.00. In other words, the introduction of experimental thresholds improves the positions of 
resonances with n > 3 relative to the ionization threshold but not those within n = 2, as 
previously mentioned. It may be appreciated in Table 9 that this finding is corroborated 
by the AE[2, 3) resulting from the threshold-energy shifting in McLaughlin & Kirby (1998). 
It is also seen in Table 9 that for the O isonuclear sequence AE{2, 3) ~ 2; for N > 2 and 
AE(2, 3) f« z + 1 for = 1 if A^(2, 3) is expressed in Ryd. 

Regarding K-resonance widths in O I, from the Auger rates listed in Table 8 widths of 
159 meV (AS2) and 163 meV (HFl) are obtained for [ls]2p^ "^P". These are in good accord 
(~15%) with the measurements of 140 meV (Krause 1994; Menzel et al. 1996) and 142 eV 
(Stolte et al. 1997) and the theoretical estimates of 169 meV (Saha 1994) and 139 meV 
(Petrini & Araujo 1994). The value of 185.22 meV hsted by McLaughhn & Kirby (1998), 
on the other hand, is somewhat high. These authors also predict for the [ls]2p^('*P)np ^P° 
widths a decrease with n. As previously mentioned, the utter dominance of KLL decay (see 
Eqs 7-9) gives rise to constant resonance widths confirmed with ASl for 3 < n < 5 to be at 
103 meV and with KLn branching ratios of 5 parts in 1000 or less. 
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7. Opacities 

In order to investigate the variations of the oxygen K-edge structure with plasma condi- 
tions, our best atomic data for O ions have been incorporated in the XSTAR modelhng code 
(Kallman & Bautista 2001). Opacities for ionization parameters in the range 0.001 < ^ < 10 
have then been obtained with a simple model denoted by the following specifications: gas 
density of 10^^ cm~^; X-ray luminosity source of 10^^ erg/s; and solar abundances for H, He, 
and O. Since the density is assumed constant, the well-known expression for the ionization 
parameter 

e = L/nR^ (12) 

is implemented where L is the luminosity of the source, R its distance, and n the density of 
the gas (Tarter et al. 1969). Therefore, the monochromatic opacities are calculated assuming 
that ionization and heating are mainly due to an external source of continuum photons. 

In Fig. 8 opacities for different values of log ^ are plotted in the 400-800 eV energy range 
showing a rich variability mainly determined by the resonance structure and K edges of the 
abundant ionic species. For a lowly ionized plasma (log^ = —3), the opacity is dominated 
by the O I edge at ~545 eV and the K« resonance structure belonging to O I and O II in 
the 530-540 eV energy interval. For log,^ = —2 the O II edge is also conspicuous at '^565 
eV, but by log,^ = —1.5 is almost vanished being replaced by the O ill and O iv edges at 
~600 eV and ~630 eV, respectively. While the latter is still distinctive for log^ = — 1, the 
O V edge arises at ~665 eV and becomes the dominant feature for log^ = 0. At a high 
degree of ionization (log^ = 1), a three-edge structure (O V, O VI, and O Vll) is observed 
together with strong absorption features due to K/3 and high-n resonances that in general 
enhance edge smearing, the opposite characteristics of low ionization. This behavior at the 
two extremes of the ionization parameter is similar to that reported for the Fe isonuclear 
sequence by Palmeri et al. (2002), an important difference being that the oxygen opacities 
show important absorption features in the edge region due to the always open L shell. 



8. Discussion and conclusions 

Extensive computations of the atomic data required for the spectral modelling of K- 
shell photoabsorption of oxygen ions have been carried out in a multi-code approach. This 
method has proven to be — as in a previous systematic treatment of the Fe isonuclear sequence 
(Bautista et al. 2003; Palmeri et al. 2003a,b; Bautista et al. 2004; Mendoza et al. 2004) — 
effective in sizing up the relevant interactions and in obtaining accuracy rating estimates, 
the latter often overlooked in theoretical work. Inter-complex CI and relativistic corrections 
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are found to be small, but CRE are shown to be conspicuous in most processes involving 
K- vacancy states. The inclusion of CRE in ion representations is therefore compelling which 
in practice is managed (e.g. in the HFR and AUTOSTRUCTURE codes) by considering non- 
orthogonal orbital bases for K-vacancy and valence states; otherwise (e.g. BPRM) orthogonal 
orbital bases are optimized by minimizing energy sums that span both types of states. 

By considering several approximations and comparing with reported experimental and 
theoretical data, the present level energies and wavelength for O ions with electron occupan- 
cies N < A are estimated to be accurate to within 0.5 eV and 0.02 A, respectively. For ions 
with > 4, the absence of measurements in particular for 5 < < 6, the wide experimental 
scatter for > 6, and discrepancies among theoretical data sets make the estimate of reli- 
able accuracy ratings less certain. It is found that ior N > A our level energies (wavelengths) 
are probably somewhat high (short) by up to 1.5 eV (0.05 A) in neutral oxygen. Our best 
results are from approximation HFl. 

The radiative data {gf -values and total widths) computed with approximations AS2 and 
HFl are expected to be accurate to 10% except in O l where the rating is not better than 
25% and in transitions involving strongly mixed levels, namely the [ls]2p^ ^^3/2 and ^P3/2 
levels in O ll and to a lesser extent the four [ls]2s2p ^Pj levels in O VI, where it is unrehable. 

Admixture in the latter levels also affect Auger rates which otherwise are expected to be 
accurate to 20%. The complete and ranked set of radiative and Auger data for transitions 
involving K-vacancy levels in O ions is one of the main achievements of the present work. 

A second contribution is the high-energy photoabsorption cross sections for ions with 

N > 3 that depict in detail the resonance structure of the K edge. Such structure has been 
shown to be dominated by KLL damping (often neglected in previous work) which gives the 
edge a blunt imprint with arguably diagnostic potential. A comparison with the experimental 
K resonance positions in O I results in conclusive statements about the importance of R- 
matrix computations with well balanced close-coupling expansions, i.e. where both terms 
of expression (4) are separately convergent. In this respect, the accuracy of the energy 
interval between the n — 2 and n — 3 components of a K resonance series gives a measure 
of the close-coupling balance, and for the O isonuclear sequence, this energy interval when 
expressed in Ryd has been shown to be close to the effective charge, z = Z — N + 1, for 
members with N >2 and near z + 1 for = 1. This finding could be exploited in spectral 
identification. 

The present atomic data sets have been incorporated in the XSTAR modelling code in 
order to generate improved opacities in the oxygen K-edge region using simple constant 
density slab models. Their behavior as a function of ionization parameter, shown in detail 
for the first time, is similar to that of the iron K edge with the exception of omnipresent 
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Ka resonances. Finally, as a continuation to the present work we intend to treat with 
the same systematic methodology the Nc isonuclear sequence and then other sequences of 
astrophysical interest (Mg, Si, S, Ar, and Ca), always bearing in mind the scanty availability 
of laboratory measurements which as shown in the present work are unreplaceable requisites 
in theoretical fine-tuning. 

MAB acknowledges partial support from FONACIT, Venezuela, under contract No. Sl- 
20011000912. TWG was supported in part by the NASA Astronomy and Physics Research 
and Analysis Program. PP was supported by a return grant of the Belgian Federal Scientific 
Policy. 
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Fig. 1. — Energy differences for K- vacancy levels in the O isonuclear sequence between those 
computed with approximation AS2 and those with ASl (open squares), HFl (crosses) and 
RMl (filled circles). Differences with ASl and RMl are due mainly to CRE. The agreement 
between AS2 and HFl is within ±0.5 eV. 
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Fig. 2. — Wavelength differences for Ka transitions in the O isoelectronic sequence between 
those computed with approximation AS2 and those with ASl (open squares), HFl (crosses), 
and RMl (sohd circles). The agreement between AS2 and HFl is within ±0.02 A. 
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Fig. 3. — Comparison of gf -values for Ka transitions in O ions computed with approxima- 
tions AS2 and HFl. It may be seen that, while the agreement for most transitions is within 
10%, the HFl g^Z-values for O I are consistently larger by 25%. 
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Fig. 4. — Comparison of Auger widths (s^^) for K-vacancy levels in O ions computed with 
approximation AS2 and those with ASl (filled circles) and HFl (open circles). Differences 
between ASl and AS2 (as large as 35%) are due to CRE. It may also be seen that HFl are 
on average 10% lower than AS2. 
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Fig. 5. — High-energy photoionization cross sections of O ions showing the structure of the K 
edge. Sohd curve — RMl; dotted curve — Pradhan et aL (2003); broken curve — Reilman 
& Manson (1979). 
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Fig. 6.— Quantum defects for the (a) [ls]2p^(^P)np and (b) [ls]2p4(2p)np 3p° (3 < n < 
6) resonance series of O I plotted relative to their respective threshold energies. Circles — 
present RMl results. Filled circles — experiment of Stolte et al. (1997). Filled triangles — 
i?-matrix data by McLaughlin & Kirby (1998). 
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Fig. 7.— Quantum defects for the (a) [ls]2p^(4p)np and (b) [ls]2p^(2p)np (3 < 
n < 10) resonance series of O I plotted relative to their respective threshold energies. Filled 
squares — experiment of Krause (1994) and Menzel et al. (1996). Filled circles — experiment 
of Stolte et al. (1997). Triangles — i?-matrix quantum defects of McLaughlin & Kirby 
(1998) relative to the experimental thresholds of Krause (1994). Filled triangles — i?-matrix 
quantum defects of McLaughhn & Kirby (1998) relative to the experimental threshold by 
Stolte et al. (1997). 
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Fig. 8. — Opacities for a photoionized gas in the oxygen K-edge region for ionization param- 
eters in the range — 3 < log^ < 1. 
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Table 1. Calculated energies (eV) for valence and K- vacancy levels in O ions 



N 


i 


Level 


ASl 


AS2 


HFl 


RMl 


N 


i 


Level 


ASl 


AS2 


HFl 


RMl 


2 
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0.000 


0.000 


0.000 
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I 


ZP Dgyj 


0.000 
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0.000 
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2 


2 


ri„io 3 no 

[lsJ2p "Fj^ 
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568.1 




7 


2 


o„3 2t^o 
2P Dg/2 


3.972 


3.711 


3.770 


3.615 


2 


3 


[is]2p ip; 
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573.5 


573.7 




7 


3 


2p3 2d°^2 


3.975 


3.714 


3.770 


3.614 


3 


1 


2s 2Si/2 


0.000 


0.000 


0.000 


0.000 


7 


4 


2p3 2po^^ 


5.393 


5.064 


5.131 


5.655 


3 


2 


[ls]2s2p 2po^2 


561.7 


563.0 


562.7 


563.4 


7 


5 


2p3 2po^^ 


5.395 


5.066 


5.131 


5.654 


3 


3 


[ls]2s2p 2p°^2 


561.7 


563.1 


562.7 


563.4 


7 


6 


[ls]2p4 4P5/2 


533.2 


531.8 


531.7 


532.5 


3 


4 


[Isl2s2p 2p?,^ 
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568.6 


568.0 


568.0 


7 


7 


[lsl2p-* 4Po/2 
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531.9 
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3 
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[±bjz&zp 3/2 
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568.1 


568.0 


7 


g 


[ibjzp ^1/2 
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4 


1 


ZS bo 


0.000 


0.000 


0.000 
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7 


9 


[Isjip D5/2 
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OOD.l 
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4 
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[ls]2p ipj 
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554.4 
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7 


10 
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5 
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2p 2po^^ 
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7 
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537.4 
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536.3 
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538.3 


540.4 


5 


4 


[ls]2p2 2D3/2 


544.1 


544.6 


544.3 


544.4 


8 


1 


2p4 3P2 


0.000 


0.000 


0.000 


0.000 


5 


5 


[ls]2p2 2p^/2 


544.8 


545.2 


545.0 


545.3 


8 


2 


2p'' 3p^ 


0.023 


0.020 


0.018 


0.000 


5 


6 


[ls]2p2 2P3/, 


544.8 


545.3 


545.0 


545.3 


8 


3 


2p4 3po 


0.033 


0.029 


0.027 


0.000 


5 


7 


[ls]2p2 2Si/2 


546.5 


547.1 


547.3 


548.2 


8 


4 


2p4 ID2 


2.353 


2.181 


2.200 


2.149 


6 


1 


2p2 3po 


0.000 


0.000 


0.000 


0.000 


8 


5 


2p4 iSo 


4.244 


3.975 


3.979 


4.833 


6 


2 


2p2 


0.016 


0.014 


0.014 


0.014 


8 


6 


[ls]2p-''' 3po 


531.4 


528.8 


528.2 


529.7 


6 


3 


2p2 3P2 


0.044 


0.037 


0.040 


0.040 


8 


7 


[ls]2p'"'' 3po 


531.4 


528.8 


528.2 


529.7 


6 


4 


2p2 ID2 


2.958 


2.750 


2.848 


2.711 


8 


8 


[ls]2p5 3po 


531.4 


528.9 


528.3 


529.7 


6 


5 


2p2 iSo 


5.390 


5.000 


5.237 


6.089 


8 


9 


[ls]2p5 ipj 


533.5 


531.2 


530.8 


532.7 


6 


6 


[ls]2p3 3d° 


537.5 


536.9 


537.0 


536.9 
















6 


7 


[ls]2p3 3d° 


537.5 


536.9 


537.0 


536.9 
















6 


8 


[ls]2p3 3 DO 


537.5 


536.9 


537.0 


537.0 
















6 


9 


[ls]2p3 355 


538.0 


537.4 


537.7 


538.0 
















6 


10 


[ls]2p3 3po 


538.9 


538.4 


538.6 


539.4 
















6 


11 


[ls]2p3 3P° 


538.9 


538.5 


538.7 


539.4 
















6 


12 


[ls]2p3 3p° 


538.9 


538.5 


538.7 


539.4 
















6 


13 


[ls]2p3 iD° 


540.8 


540.3 


540.6 


540.7 
















6 


14 


[ls]2p3 ipj 


542.2 
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Table 2. Comparison of valence level energies (eV) 



N Level Expt'' ASl AS2 HFl RMl 



2P 'Pv2 

2p2 

2p2 

2p2 ID2 
2p2 iSo 
2p3 4s°/^ 
2p-^ 2^0 



/2 



2p'^ ^DJ 
2p3 2po^^ 



2p 



3/2 
3 2po 



1/2 



2p4 
2p4 
2p4 

2p4 ID2 
2p4 iSo 



0.000 
0.048 
0.000 
0.014 
0.038 
2.514 
5.354 
0.000 
3.324 
3.327 
5.017 
5.018 
0.000 
0.020 
0.028 
1.967 
4.190 



0.000 
0.057 
0.000 
0.016 
0.044 
2.958 
5.390 
0.000 
3.972 
3.975 
5.393 
5.395 
0.000 
0.023 
0.033 
2.353 
4.244 



0.000 
0.049 
0.000 
0.014 
0.037 
2.750 
5.000 
0.000 
3.711 
3.714 
5.064 
5.066 
0.000 
0.020 
0.029 
2.181 
3.975 



0.000 

0.047 
0.000 
0.014 
0.040 
2.848 
5.237 
0.000 
3.770 
3.770 
5.131 
5.131 
0.000 
0.018 
0.027 
2.200 
3.979 



0.000 
0.047 
0.000 
0.014 
0.040 
2.711 
6.089 
0.000 
3.615 
3.614 
5.655 
5.654 
0.000 
0.000 
0.000 
2.149 
4.833 



■Spectroscopic tables (Moore 1998) 
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Tablc 3. Comparison of K- vacancy level energies* (eV) 



N 


State 


AS2 


HFl 




Expt 


Theory 


2 


[ls]2p 


567.7 


568.1 


568.6'' 




568.2" 


2 


[ls]2p 1P° 


573.5 


573.7 


574.0^ 




573.6" 


3 


[Is]2s2p(-^P°) 2p° 


563.1 


562.7 


562.6=^ 




562.3' 


3 


[ls]2s2p(ip°) 2po 


568.6 


568.1 


568.2'>- 




567.0' 


4 


[ls]2p ipj 


554.5 


554.4 






554.6', 553.15s 


5 


[ls]2p2 2D3/2 


544.6 


544.3 






545.8' 


5 


[ls]2p2 2P 


545.3 


545.0 






547.0' 


5 


[ls]2p2 2Si/2 


547.1 


547.3 






552.1' 


6 


[ls]2p3 305- 


536.9 


537.0 






537.2' 


6 


[ls]2p3 3s° 


537.4 


537.7 






538.6' 


6 


[ls]2p3 3pj 


538.5 


538.7 






540.7' 


7 


[ls]2p4 4p 


531.9 


531.7 


530.8(3)'', 


530.41(4)"= 


532.9', 531.0'', 534.1' 


7 


[ls]2p4 2d 


536.1 


536.1 






535.3^, 538.7' 


7 


[ls]2p4 2p 


536.2 


536.3 


535.7(3)'', 


535.23(4)-= 


535.9'', 539.2' 


7 


[ls]2p4 2 s 


538.4 


538.3 






538.1'', 542.9' 


8 


[ls]2p5 3po 


528.8 


528.2 


527.2(3)'', 


526.79(4)=, 527.85(10)'' 


528.8', 528.22', 528.33' 



* Relative to each ion ground level 

"Spectroscopic tables (Moore 1998) 

''Auger electron spectrometry (Krause 1994; Caldwell et al. 1994) 
'=Photoionization experiment (Stolte et al. 1997) 
''Photoabsorption measurements (Menzel et al. 1996) 
"Non-relativistic CI calculation (Charro et al. 2000) 
'i?- matrix calculation by Pradhan et al. (2003) 
SMulticonfiguration Dirac-Fock calculation by Chen (1985) 
''il-matrix calculation (Gorczyca & McLaughlin 2000) 
'il-matrix calculation (McLaughlin &c Kirby 1998) 
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Table 4. Calculated wavelengths and gf -values for Ka transitions 











AS2 




HFl 










AS2 




HFl 


N 


i 


J 


A 

(A) 


if J 13 


A 

(A) 


UJlJ 


N 


i 


J 


A 

(A) 




A 

(A) 


O fi T 


2 


1 


2 


21.84 


1.22E-4 


21.82 


7.35E-5 


6 


4 


13 


23.06 


1.35E+0 


23.06 


1.35E+0 


2 


1 


3 


21.62 


8.27E-1 


21.61 


7.96E-1 


6 


4 


14 


23.00 


5.09E-1 


22.99 


4.98E-1 


3 


1 


2 


22.02 


4.93E-1 


22.04 


4.37E-1 


6 


5 


14 


23.09 


3.37E-1 


23.09 


3.41E-1 


3 


1 


3 


22.02 


9.95E-1 


22.03 


8.81E-1 


7 


1 


6 


23.31 


4.21E-1 


23.32 


4.56E-1 


3 


1 


4 


21.80 


4.53E-2 


21.83 


4.40E-2 


7 


1 


7 


23.31 


2.83E-1 


23.32 


3.04E-1 


3 


1 


5 


21.80 


8.23E-2 


21.83 


8.09E-2 


7 


1 


8 


23.31 


1.41E-1 


23.32 


1.52E-1 


4 


1 


2 


22.36 


7.08E-1 


22.37 


6.52E-1 


7 


2 


9 


23.29 


6.22E-1 


23.29 


6.70E-1 


5 


1 


4 


22.77 


3.53E-1 


22.78 


3.33E-1 


7 


2 


10 


23.29 


l.OOE-1 


23.29 


1.45E-2 


5 


1 


5 


22.74 


4.36E-1 


22.75 


4.09E-1 


7 


2 


11 


23.29 


6.72E-1 


23.28 


8.18E-1 


5 


1 


6 


22.74 


1.89E-1 


22.75 


1.80E-1 


7 


3 


9 


23.29 


4.09E-2 


23.29 


4.43E-2 


5 


1 


7 


22.66 


7.15E-2 


22.65 


6.67E-2 


7 


3 


10 


23.29 


4.76E-1 


23.29 


4.74E-1 


5 


2 


3 


22.77 


5.90E-1 


22.78 


5.61E-1 


7 


3 


11 


23.29 


1.48E-2 


23.28 


5.42E-2 


5 


2 


4 


22.77 


4.21E-2 


22.78 


4.24E-2 


7 


3 


12 


23.29 


3.98E-1 


23.28 


4.30E-1 


5 


2 


5 


22.74 


2.06E-1 


22.75 


1.95E-1 


7 


4 


9 


23.35 


2.00E-1 


23.35 


2.19E-1 


5 


2 


6 


22.74 


1.09E+0 


22.75 


1.0.3E+0 


7 


1 


10 


23.35 


2.04E-1 


23.35 


4.41E-2 


5 


2 


7 


22.66 


1.68E-1 


22.66 


1.52E-1 


7 


4 


11 


23.35 


1.66E-1 


23.34 


3.61E-1 


6 


1 


8 


23.09 


1.41E-1 


23.09 


1.42E-1 


7 


4 


12 


23.34 


7.48E-2 


23.34 


8.09E-2 


6 


1 


9 


23.07 


1.30E-1 


23.06 


1.30E-1 


7 


4 


13 


23.25 


2.12E-1 


23.25 


2.29E-1 


6 


1 


11 


23.03 


8.40E-2 


23.02 


8.18E-2 


7 


5 


10 


23.35 


1.38E-2 


23.35 


9.82E-2 


6 


2 


7 


23.09 


3.17E-1 


23.09 


3.18E-1 


7 


5 


11 


23.35 


1.67E-1 


23.34 


l.OOE-1 


6 


2 


8 


23.09 


l.OOE-1 


23.09 


9.95E-2 


7 


5 


12 


23.34 


1.46E-1 


23.34 


1.61E-1 


6 


2 


9 


23.07 


3.82E-1 


23.06 


3.82E-1 


7 


5 


13 


23.25 


l.OlE-1 


23.25 


1.08E-1 


6 


2 


10 


23.03 


1.07E-1 


23.02 


1.05E-1 


8 


1 


6 


23.45 


4.31E-1 


23.47 


5.25E-1 


6 


2 


11 


23.03 


6.52E-2 


23.02 


6.31E-2 


8 


1 


7 


23.44 


1.44E-1 


23.47 


1.75E-1 


6 


2 


12 


23.03 


9.31E-2 


23.02 


9.12E-2 


8 


2 


6 


23.45 


1.44E-1 


23.47 


1.80E-1 


6 


3 


6 


23.10 


5.75E-1 


23.09 


5.77E-1 


8 


2 


7 


23.45 


8.61E-2 


23.47 


1.08E-1 


6 


3 


7 


23.09 


9.36E-2 


23.09 


9.44E-2 


8 


2 


8 


23.44 


1.15E-1 


23.47 


1.44E-1 


6 


3 


8 


23.09 


5.77E-3 


23.09 


6.05E-3 


8 


3 


7 


23.45 


1.15E-1 


23.47 


1.44E-1 


6 


3 


9 


23.07 


6.11E-1 


23.06 


6.08E-1 


8 


4 


9 


23.44 


6.28E-1 


23.45 


7.85E-1 


6 


3 


10 


23.03 


3..58E-1 


23.02 


3.51E-1 


8 


5 


9 


23.52 


1.05E-1 


23.53 


1.32E-1 


6 




li 


2:?.():-! 


l.:-!lE-i 


2:-!, 02 


1.2I)E-1 
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Table 5. Wavelength (A) comparison for Ka transitions 



N 


i 




AS2 


HFl 


Expt 


Other theory 


2 


ls2 ISo 


[ls]2p 3p° 


21.84 


21.82 


21.807^ 


21.84', 21.841 




ls2 iSo 


[ls]2p ipj 


21.62 


21.61 


21.6020(3)'' 


21.60' 


3 


2s 2s 


[ls]2s2p(3p°) 2po 


22.02 


22.03 


22.038*, 22.01(1)'=, 22.0194(16)'' 


22.02', 22.06j, 22.05'', 22.00' 




2s 


[ls]2s2p(ip°) 2p° 


21.80 


21.83 


21.82* 


21. 84j ,21.87'', 21.79' 


4 


2s2 is 


[ls]2p ip° 


22.36 


22.37 


22.38(1)=, 22.374(3)d 


22.38', 22.4lj, 22.35''^, 22.33' 


5 




[ls]2p2 2D3/2 


22.77 


22.78 


22.74(2)"= 


22.73'', 22.73' 




2p ^P°/, 


[ls]2p2 2p 


22.74 


22.75 




22.67'', 22.78' 






[ls]2p2 2Si/2 


22.66 


22.65 




22.46'', 22.73' 


6 


2p2 3Po 


[ls]2p3 3d° 


23.09 


23.09 


23.17(1)<= 


23.08'', 23.11' 




2p2 3Po 


[ls]2p3 330 


23.07 


23.06 


23.00(2)"= 


23.02'', 23.05' 




2p2 3po 


[ls]2p3 3po 


23.03 


23.02 




22.93'', 22.98' 


7 


2p3 *S° 


[ls]2p4 4p 


23.31 


23.32 


23.36(1)° 


23.27"', 23.30' 




2p3 2d° 


[ls]2p4 2p 


23.29 


23.28 


23.29(1)° 






2p3 2po 


[ls]2p4 2p 


23.35 


23.34 


23.36(1)" 




8 


2p'4 3p 


[ls]2p5 3po 


23.45 


23.47 


23.52(3)«, 23.489(5)*, 23.536(2)8 


23.45'' 



23.508(3)'' 



"Spectroscopic tables (Moore 1998) 

''Spectroscopic measurements (Engstrom & Litzen 1995) 
"^Spectroscopy of NGC 5548 (Steenbrugge et al. 2003) 

"'Electron beam ion trap measurements (Schmidt et al. 2004) 
"= Auger electron spectrometry (Krausc 1994; Caldwell et al. 1994) 
' Photoabsorption measurements (Menzel et al. 1996) 
sPhotoionization experiment (Stolte et al. 1997) 
•"ISM observations (Juett et al. 2004) 

'1/Z expansion method (Vainshtein & Safronova 1971, 1978) 
JMulticonfiguration Dirac-Fock calculation (Chen 1985, 1986) 
''iZ-matrix calculation (Pradhan et al. 2003) 
'hullac calculation (Behar & Kahn 2002) 
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Table 6. Comparison of theoretical /-values for Ka transitions 



N 


i 


j 


AS2 


HFl 


RM2'' 


3 


2s 2S 


[ls]2s2p(3p°) 2po 


0.744 


0.659 


0.576 


3 


2s 2S 


[ls]2s2p(lp°) 2po 


0.077 


0.062 


0.061 


4 


2s2 is 


[ls]2p lp° 


0.708 


0.652 


0.565 


5 




[ls]2p2 2D3/2 


0.177 


0.167 


0.132 


5 




[ls]2p2 2p 


0.313 


0.295 


0.252 


5 


2P 'P?/2 


[ls]2p2 2Si/2 


0.036 


0.033 


0.027 


6 


2p2 


[ls]2p3 


0.141 


0.142 


0.119 


6 


2p2 3Po 


[ls]2p3 3S° 


0.130 


0.130 


0.102 


6 


2p2 


[ls]2p3 3po 


0.084 


0.082 


0.067 


7 


2p3 4Sc^^ 


[ls]2p4 4p 


0.212 


0.228 


0.184 


8 


2p4 3P2 


[ls]2p5 3p 


0.115 


0.140 


0.113 



''ii-matrix data by Pradhan et al. (2003) 



Table 7. Comparison of branching ratios for KLL Auger transitions 



N 


i 


j 


ASl 


AS2 


HFl 


Expf^ 


MCDF'' 




7 


[ls]2p4 4p 




0.42 


0.42 


0.43 


0.55(13) 


0.36 


0.43 






[2s]2p3 


0.41 


0.40 


0.40 


0..33(3) 


0.46 


0.42 






[2s2]2p4 


0.17 


0.18 


0.18 


0.12(2) 


0.18 


0.15 


7 


[ls]2p4 2p 


2p2 


0.53 


0.53 


0.51 


0.53(9) 


0.45 


0.47 






[2s]2p3 


0.27 


0.27 


0.28 


0.28(4) 


0.31 


0.28 






[2s2]2p4 


0.20 


0.20 


0.21 


0.19(3) 


0.24 


0.25 


8 


[ls]2p5 3po 


2p3 


0.53 


0.51 


0.55 


0.606 










[2s] 2?" 


0.35 


0.35 


0.32 


0.30 










[2s2]2p5 


0.12 


0.14 


0.13 


0.094 







"Auger electron spectrometry (Caldwell & Krause 1993; Caldwell et al. 1994) 
t'MCDF calculation by M. H. Chen as quoted in Caldwell et al. (1994) 
"^Close-coupling calculation (Petrini & Araiijo 1994) 



Table 8. Calculated radiative and Auger widths for K-vacancy levels 



AS2 HFl 



N 


i 


Level 




















(s-1) 


(s-1) 


2 


2 


[ls]2p 3p° 


6.60E+08 




3.44E+08 




2 


3 


[ls]2p ip° 


3.93E+12 




3.80E+12 




3 


2 


[ls]2s2p ^P°_/2 


3.39E+12 


6.60E+12 


3.00E+12 


7.53E+12 


3 


3 


[ls]2s2p 2p°^2 


3.42E+12 


5.76E+12 


3.02E+12 


6.96E+12 


3 


4 


[ls]2s2p 2p°^2 


3.19E+11 


7.56E+13 


3.12E+11 


6.67E+13 


3 


5 


[ls]2s2p 2po^2 


2.90E+11 


7.60E+13 


2.87E+11 


6.73E+13 


4 


2 


flsl2p ip° 

L J 1 


3.40E+12 


9.10E+13 


3.10E+12 


8.15E+13 


5 


3 


[ls]2p2 2D5/2 


1.34E+12 


2.19E+14 


1.22E+12 


1.93E+14 


5 


4 


[ls]2p2 2D3/2 


1.35E+12 


2.18E+14 


1.28E+12 


1.92E+14 


5 


5 


[ls]2p2 2p.^^2 


4.31E+12 


9.29E+13 


4.05E+12 


7.57E+13 


5 


6 


[ls]2p2 2p3^2 


4.31E+12 


9.30E+13 


4.05E+12 


7.58E+13 


5 


7 


[ls]2p2 2Si/2 


1.63E+12 


1.89E+14 


1.50E+12 


1.66E+14 


6 


6 


[ls]2p3 


1.05E+12 


2.37E+14 


1.06E+12 


2.11E+14 


6 


7 


[ls]2p3 


1.05E+12 


2.37E+14 


1.06E+12 


2.10E+14 


6 


8 


[ls]2p3 ^DJ 


1.05E+12 


2.37E+14 


1.06E+12 


2.10E+14 


6 


9 


[ls]2p3 3s° 


4.81E+12 


8.87E+13 


4.81E+12 


7.24E+13 


6 


10 


[ls]2p3 3po 


1.20E+12 


2.22E+14 


1.18E+12 


1.98E+14 


6 


11 


[ls]2p3 3po 


1.20E+12 


2.22E+14 


1.18E+12 


1.98E+14 


6 


12 


[ls]2p3 3po 


1.20E+12 


2.22E+14 


1.18E+12 


1.98E+14 


6 


13 


[ls]2p3 lD§ 


3.46E+12 


2.05E+14 


3.46E+12 


1.84E+14 


6 


14 


[ls]2p3 ipj 


3.62E+12 


1.90E+14 


3.60E+12 


1.71E+14 


7 


6 


[ls]2p4 4p,./2 


8.67E+11 


2.25E+14 


9.35E+11 


2.04E+14 


7 


7 


[ls]2p4 4p3/2 


8.67E+11 


2.25E+14 


9.35E+11 


2.04E+14 


7 


8 


[ls]2p4 4pj/2 


8.68E+11 


2.25E+14 


9.35E+11 


2.04E+14 


7 


9 


[ls]2p4 2D5/2 


1.78E+12 


2.60E+14 


1.93E+12 


2.44E+14 


7 


10 


[ls]2p4 2D3/2 


2.46E+12 


2.36E+14 


1.96E+12 


2.43E+14 


7 


11 


[ls]2p4 2P3/2 


3.17E+12 


2.10E+14 


4.17E+12 


1.73E+14 


7 


12 


[ls]2p4 2p^/2 


3.85E+12 


1.86E+14 


4.20E+12 


1.72E+14 


7 


13 


[ls]2p4 2Si/2 


1.96E+12 


2.43E+14 


2.10E+12 


2.28E+14 


8 


6 


[ls]2p-"'' 3po 


1.39E+12 


2.41E+14 


1.70E+12 


2.48E+14 


8 


7 


[ls]2p5 3po 


1.39E+12 


2.41E+14 


1.70E+12 


2.48E+14 


8 


8 


[ls]2p5 3P° 


1.39E+12 


2.41E+14 


1.70E+12 


2.48E+14 


8 


9 


[ls]2p5 lp° 


2.99E+12 


2.20E+14 


3.66E+12 


2.29E+14 
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Table 9. Comparison of A£^(2, 3)* for O ions 



N 


z 


Present 


Expt 


Other theory 


8 


1 


1.06 


1.06'', 0.99"= 


1.02«, 1.00* 


7 


2 


1.93 






6 


3 


2.90 






5 


4 


3.73 






4 


5 


5.03 






3 


6 


5.64 






2 


7 


7.00 


7.02^ 




1 


8 


8.90 


8.90<* 





* Energy interval in Ryd between the lowest n = 2 
and n = 3 resonances 

"Approximation RMl for Af > 2 and ASl for TV < 2 

''Photoionization experiment (Stolte et al. 1997) 

"^Photoabsorption measurements (Menzel et al. 
1996) 

''Spectroscopic tables (Moore 1998) 

"^R-matrix calculation (McLaughlin & Kirby 1998) 
with the experimental thresholds of Krause (1994) 

'R-matrix calculation (McLaughlin & Kirby 1998) 
with the experimental thresholds of Stolte et al. 
(1997) 



